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We study primordial fluctuations generated during inflation in a class of models motivated by 
the DBI Galileons, which are extensions of the DBI action that yield second order field equations. 
This class of models generalises the DBI Galileons in a similar way with K-inflation. We calculate 
the primordial non-Gaussianity from the bispectrum of the curvature perturbations at leading order 
in the slow-varying approximations. We show that the estimator for the equilateral-type non- 
Gaussianity, /nl"*' i can be applied to measure the amplitude of the primordial bispectrum even in 
the presence of the Galileon-like term although it gives a slightly different momentum dependence 
from K-inflation models. For the DBI Galileons, we find -0.32/c 2 < /X" < -0.16/c? and large 
primordial non-Gaussianities can be obtained when c 3 is much smaller than 1 as in the usual DBI 
inflation. In G-inflation models, where a de Sitter solution is obtained without any potentials, the 
non-linear parameter is given by J^l"' = 4.62r _2,/3 where r is the tensor to scalar ratio, giving a 
stringent constraint on the model. 



I. INTRODUCTION 

The DBI inflation model [l| is one of the most interesting possibilities to realise large non-Gaussianity of the 
Cosmic Microwave Background (CMB) temperature fluctuations. Non-Gaussianity of the curvature perturbation in 
DBI inflation has been studied extensively (see also [23|, HiJ for reviews) . 

Recently, a very interesting extension of the DBI inflation model, so-called "DBI Galileons" was proposed by de 
Rham and Tolley [25| (see also [HI). This is based on the relativistic extension of the Galileon model [27l - l29j (for 
studies of cosmology based on the Galileon field, see Refs. [3fil - |40l |). The simplest example is a single field model 
that arises from a probe brane action in the five-dimensional spacetime. Let's consider the following four-dimensional 
induced action on the probe brane 

S = J d^x^—g (A - MlK) , (1) 

where A is a tension of the brane, Mg is the five dimensional Planck constant, g^u is the induced metric on the brane, 

9nu = + dvird^ir, (2) 

and K is a trace of the extrinsic curvature K u 



VI + (<9tt) 2 

Here it is a modulus describing the position of the bra ne. Using the fact that the inverse metric is given by g^ v = 
if" — 7 2 <9 M 7r9 ly 7r where 7 is the Lorentz factor 7 = l/\/l + (dir) 2 , the brane action is written as 

S = -A J d^xyjl + (<9tt) 2 + Ml J d 4 x (Dtt - 7 2 d A A7rc^7n9 ,y 7r) . (4) 

By integrating by part and discarding the total derivative terms, this action can be rewritten as 

S = - X J d 4 x^/TTWy + J d 4 x (j 2 (dn) 2 nTr + cy 7 2 )(c>V)(c)7r) 2 ) . (5) 

The first term is the usual DBI action. The higher order terms look containing the higher derivatives but the equation 
of motion is at most second order in derivatives. In the non- relativistic limit 7 — > 1, this reduces the Galileon model 
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where the action is invariant under the Galileon symmetry 9 M 7r — ¥ d^n + c^. It is also possible to include two more 
higher order terms in the action but in this paper we focus our attention to the leading order cubic order terms. 

Recently, Refs. [4l[ and [13] considered a generalisation of this model. This generalisation is based on the extension 
of DBI inflation to K-inflation. The generalised action is given by 

S = J d A x^—g{P{X, 4>) - G(X, <t>)U<t>) , (6) 

where (j> is the same degree of the freedom as 7T, but it is defined so that (j> has a dimension of mass and 
X = —(l/2)(dcf)) 2 . In Eq. ([6]), P(cj),X) and G(<p,X) are arbitrary functions of (f> and X. Precisely speaking this 
generalisation does not include the action ([5]). However, in the DBI inflation, the Lorentz factor, 7, varies very slowly 
and at leading order in the slow-varying parameters, which is usually used to calculate the leading order contribution 
to non-Gaussianity, the last term in ([5]) does not play a role. Thus the action ([6|) is general enough to include the 
case of the DBI Galileons. Again the equation of motion is at most second order in derivatives. 

In this paper, we study primordial fluctuations generated during inflation described by the action ©. We calculate 
the power spectrum of the curvature perturbation as well as the bispectrum of the curvature perturbations. Two 
examples of the model are considered: one is the DBI Galileons described by the action ([5|). The other is G-inflation 
models proposed by Ref. [42|. This model is based on a specific choice of the functions P(X) and G(X) that realises 
a de Sitter solution without any potentials. 

This paper is organised as follows. In section II, we introduce a model studied in this paper. The power spectrum 
is calculated in section II and the bispectrum is discussed in section III. In section IV, we consider two examples, the 
DBI Galileons and G-inflation. Section V is devoted to conclusions. 



II. MODEL 



Assuming that (j> described by the action © is minimally coupled to gravity, we consider a class of models described 
by the following action: 

S = \ J d 4 xV=d{M 2 pl R + 2P{<P, X) - 2G(0, X)D(f>] , (7) 

where M p i is the Planck mass. In the background, we are interested in flat, homogeneous and isotropic Friedmann- 
Robcrtson- Walker universes described by the line element 

ds 2 = -dt 2 + a 2 {t)8 l] dx t dx j , (8) 

where a(t) is the scale factor. It can be shown that the energy density and pressure of the field are given by 

p = 2P X X - P + 6G X H j>X - 2G^X , (9) 
V = P-2(G, <t> + G, x 4>)X, (10) 

where H — a /a is the Hubble parameter, a dot represents a derivative with respect to cosmic time t and the subscripts 
t x and y rj, denote derivatives with respect to X and <f>, respectively. The Friedmann equation and the field equation 
are given by 

ZM 2 pl H 2 =p, (11) 
P x (4> + 3#</>) + 2P XX X4> + 2P X *X - P A - 2G, (0 + 3ff0) - 2G^X{if> - 3ff0) 
+6G tX [(HXy + 3H 2 X] - 2G^ + 6G, XX HXX = . (12) 

It is useful to define a slow- varying parameter 

H XP X + ZG, X H4>X 



H 2 M 2 H 2 



(13) 



where for the second equality, we have assumed the quantities \4>/(H<j))\ and \G^cf)/(GH)\ are much smaller than 1. 
Since we are interested in fluctuations generated during inflation, we will consider the background that satisfies the 
slow-varying conditions which are given by |e| <§; 1, together with \(j)/(H(j))\ <C 1 and \G^(j)/(GH)\ <C 1. 



3 



III. POWER SPECTRUM 



We are interested in the primordial curvature perturbation on uniform density hypersurfaces, £, on large scales, 
which is directly related to temperature anisotropies in the Cosmic Microwave Background (CMB). In order to 
calculate the statistical quantities of £ at leading order in the slow-varying approximations, we first calculate the 
bispectrum of the fluctuations of inflaton (j) in the flat gauge where the three-dimensional metric takes the form 
hij = a 2 5ij, and then relate it to that of £ using the relation obtained from the delta-N formalism (43l. Eij 



(14) 



In Eq. (p~4|) . (f> is the background value and Q is the perturbation in the flat gauge. In this gauge, at leading order in 
the slow- varying approximations, the second-order action is expressed as 



.s'. - J dtd 3 xa 3 PxSXW + X _ Pxx ^SX (1 ^ 2 - G, X 8X^U^ - G, x 6X^n^ - X -G, xx [bX^ n<j>™ 



with 



5X™ = i SX^ = X -Q 2 - ^3 l Qd t Q , □0 (o) = -'dH4> , D^ (1) = -Q - 3HQ + ^d l d v 
Introducing the sound speed for the scalar perturbations 

2 P x + 44>HG,x 



Px + 2XP X x + QH<t>{G, x + XG, XX ) 
and integrating by parts, we can write the second order action as 



S 2 = j dtd A x^(P,x+±4>HG, x ) 



q 2 - ^d l Qd % 



The perturbations in the interaction picture are promoted to the quantum operators as 

Q(r,x) = -i TT / d 3 kQ(r, k)e ak x , Q(r, k) = u(t, k)a(k) + u*(t, -k) t(-k) , 

(2ny J 



,(15) 



(16) 



(17) 



(18) 



(19) 



where a(k) and a^(— k) are the annihilation and creation operator respectively. They satisfy the usual commutation 
relations 



[a(k x ),at(k 2 )] =(2^) 3 ,5( 3 )(k 1 -k 2 ), [ a (k 1 ),a(k a )] - [ot(k 1 ),ot(k a )] =0. 

From the second order action (|18p . the solution for the mode functions is given by 

H 1 



u(t, k) = 



■ X ) 

It is convenient to introduce the following parameters 

G^X 
V ~ M* t H 



k 3/: 



(1 + ikc s r) e 



(20) 



(21) 



e = e + v . 



(22) 



where e coincides with e when there is no Galileon-like term, e is also much smaller than 1 for ^ « 1. Then, the 
power spectrum of Q and £ are given by 



(Q(k 1 )Q(k 2 )) = (2 7 r) 3 <5( 3 )(k 1 + k 2 )P Q ^ , V Q = 47r2 * 2 g . 



2tt 2 



(c(koc(k 2 )) = (2^) 3 ^ 3 )(k! + k 2 )^^r . n 



1 H 2 



(23) 
(24) 
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which are evaluated at the time of the sound horizon exit, c s k = aH. It is worth noting that it is e not e which appears 
in Eq. (fM|) . which gives the behaviour of the power spectrum that is different from the usual K-inflation model as we 
will see below. Defining additional slow- varying parameters 

^w- (25) 

the spectral index of the primordial power spectrum is given by 

n s -l= C 7 V ' = -2e-f)-s. (26) 

We need to require e, 77, s to be very small in order to realise the almost scale invariant power spectrum. We have 
confirmed that this result is consistent with the one obtained in Ref. [42| when the conditions \(f>/(H(f>)\ <C 1 and 
\G^(j)/(GH) \ <C 1 are satisfied. Notice that fj is different from the usual 77 defined by 7/ = e/(eH). 
The power spectrum and spectral index of tensor perturbations are given by the usual expression 

VT = ^ nT = ~ 2e - (27) 

In the usual K-inflation, and tensor to scalar ratio r = Vt/'Pc, are n °t independent, and there is a so-called 
"consistency relation" r — ~%c s tit [45(. However, it is clear that this relation does not hold in the presence of the 
Galileon-like term. Instead, we have 

r = -8c s {n T - 2v) , (28) 

where v is given by Eq. 



IV. BISPECTRUM 

The third order action can be obtained in the same way as 

To n. . . . 1 

(29) 



n 3 

So = I dtd 3 x^- 



C\Q 3 + ^Q&QdiQ + -%d l Qd l Q&d J Q + -^Qd i Qd l Q 



where 



C x = ^X 2 P XXX + XP XX + 2H<f>X 2 G, xxx + 5Hcj>XG, xx + H^G, X , 

C 2 = - (XP XX + 3H4>XG^ XX + Hd>G, x ) , C* 3 = ^Hj>G, x , C 4 = 2Hj>G, x + 2H^XG, XX . (30) 
The vacuum expectation value of the three point operator in the interaction picture is written as [43l . |46| 

(Q(«,k 1 )Q(i,k 2 )Q(i ) k 3 )) = -i f <fi'([Q(t ) k 1 )Q(i,k 2 )Q(t,kB),H/(*)]>, (31) 

Jto 

where to is some early time during inflation when the field's vacuum fluctuation are deep inside the sound horizon 
and t is some time after the sound horizon exit. If one uses a conformal time, it is a good approximation to perform 
the integration from —00 to because r rj —(aH)^ 1 . Hj denotes the interaction Hamiltonian and it is given by 
Hi = — L3, where L3 is the lagrangian obtained from the action (|29p . Using the solution for the mode function and 
the commutation relations for the creation and annihilation operators, we get 

(Q(k 1 )Q(k 2 )Q(k 3 )> = (2 7 r) 3 ^ 3 )(k 1 + k 2 + k 3 ) — fl v^na 1 ^ > ( 32 ) 

(Pa- + 44>HG tX ) d (p *-h=\ K i 



where 



= 3 I Cx - % )Ai + ^A 2 + §U 3 , (33) 
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and Ci, C 2 , C3, C4 and C5 are given by Eq. (l30l) . In Eq. (|33|) . we have introduced the shape functions Ai, A2 and 
A3 as 

1.21.21,2 

A - , (34) 

A 2 = (l + JL_H + + 2 perms. , (35) 

A 3 = %^ (l + fc ifr + ffi + *3*i + 3 ^ + 2 permg (36) 

where K = ki + k 2 + k 3 . The shapes _4i and ^2 appear in the usual K-inflation models Q and their amplitudes can 
be measured by the estimator for the equilateral-type non-Gaussianity, f^ il [47j - On the other hand, the shape ^3 
is completely new that arises from the Galileon-like term. 

For the bispectrum of f, making use of Eqs. (|14l) and ([2~4")l . we obtain the following expression: 



where 



of using the estimator for equilateral type non-Gaussianity, f^T" , for the bispectrum even in the presence of the new 



(C(k 1 )C(k 2 )C(k 3 )> = (2^) 3 < 5< 3 )(k 1 + k 2 + k 3 )(V c ) 2 F(ki, fc 2 , fc 3 ) , (37) 
%.fe,fe) = |^^/ ( ^, (38) 

f(l) = " 3 ( C 'l C s ~ C 4) ,( 2 ) = C 2 (3) = C3 

(Fx + HHG, X ) ' 2(P X + 40HG, X ) ' ' (P x + HHG, x )c 2 s 

In the following, we will study the momentum dependence of the bispectrum. Especially, we will check the validity 
f using the estimator for equilateral ty 
shape .A3. The estimator is defined by 

F^(kiMM) = (2tt) 4 (y^) (-i - i - j^ff - 737^2 + ™ + (5 P erm S .)) , (40) 

\ / \ 12 1 3 2 3 12 3 1 2 3 / 

where the permutations act only on the last term in parentheses. This shape is factorisable and it is possible to 
construct a fast optimal estimator that can be applied to the CMB map. For this purpose, it is useful to define shape 
functions F^ 1 ' (ki, k 2 , k 3 ), i = 1,2, 3 corresponding to the shapes Ai in Eq. (|38p . As mentioned before, the shape of 
the bispectrum in K-inflation is characterised by the sum of F^(h,k2,h) and F( 2 )(fc 1; k 2 , k 3 ) with /« and 
depending on P{X). For example, in DBI inflation, the relation /' 2 V/ = —2/3 holds. However, since functions 
F^(ki,k 2 ,k 3 ) and F^(ki,k 2 ,k 3 ) are not factorisable, F equll (ki, k 2 , k 3 ) is usually used to approximate the shape 
functions F^(k u k 2 ,k 3 ) and F( 2 )(fci, fc 2 , /c 3 ). 

In Figure [2 we compare F equll (k u k 2 , k 3 ) with F^Oi, fc 2 , fc 3 ), F( 2 )(fci, fc 2 , fc 3 ) and F^ih, k 2 , k 3 ) with appro- 
priate normalisations. From this, we see that not only F^^ki, k 2 , k 3 ) and F^ (ki, k 2 , k 3 ), but also the function 
F^(k 1: k 2 ,k 3 ) has a very similar shape with F eqml (ki,k 2 ,k 3 ) and it can be expected that /^j™ provides a good 
measure of the bispectrum even in the presence of the new shape function F^ 3 K We can prove this quantitatively by 
the following shape correlator C for two different shapes characterised by F and F' introduced by Ref. [48| (See also 
©113), 

C(F,F') = -= - F ( f ' F/ ) ( 41 ) 

where the overlap function J 7 is given by 

F{F,F') = j dV k F(k 1 ,k 2 ,k 3 )F'(k 1 ,k 2 ,k 3 )IlUiktwB{ki,k2,k3). (42) 

In Eq. (|4"2")l the integration is performed for the region where the triangle condition for (ki,k 2 ,k 3 ) holds and weight 
function wb is given by 

w B = -, r T ■ (43) 

" k l +k 2 + k 3 K ' 

Table U shows that the shapes F equll (k u k 2 ,k 3 ) and F^ (k u k 2 , k 3 ) are almost completely anti-correlated, which 
means that it is very reasonable to adopt the estimator to measure the amplitude of the bispectrum even in 

the presence of the shape F^(ki, k 2 ,k 3 ). 
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FIG. 1: In this group of figures, we plot the shape functions F(l,k2/ki,k 3 /ki)(k2/ki) 2 (k 3 /ki) 2 as functions of (fa/fci, k 3 /ki). 
The figures are normalised to have values f for equilateral configurations k2jk\ = k 3 /k\ = I and set to zero outside the region 
l-fe/fci < < fca/fci. We plot F (1) (fci,ifc2,ifc3), F {2) , k 2 , k 3 ) , F (3) (k 1 ,k 2 ,k 3 ) and F eg ""(fci, fc 2 , fcs) for upper left, upper 

right, lower left, lower right, respectively. 





F (1, (fa,fa,fe) 


F^ 2 \k 1 ,k 2 ,k 3 ) 


F^(k 1 ,k 2 ,k s ) 


F DBI (k 1 ,k 2 ,ks) 


Overlap 


0.936 


-0.995 


-0.99989 


-0.993 



TABLE I: Shape correlations between the factorisable equilateral shape F equ ' 1 (ki, k 2 , k 3 ) and the shapes of primordial bispec- 
tra characterised by functions F^'{ki, fe, fe), ^ &2, £3) and F^ (k\,k2,k 3 ). For comparison, we have also shown the 

correlation between the equilateral shape and that obtained in DBI inflation which is given by F i - 1 \ki,k2,k 3 ) -\- F^ 2 \ki,k2,k 3 ) 
satisfying / (2) // (1) = -2/3. 



Now that the validity of using the estimator for the equilateral-type non-Gaussianity, f^™ 1 for the bispectrum 
corresponding to the new shape A3 is confirmed, we will obtain for the bispectrum given by Eq. (|38|) . We 

match the two different shapes and ((301) so that these two shapes have the same amplitudes at the equilateral 
configuration k% = &2 = &3- This gives 

reqml = , (1) _ 85 (2) _ 65 (3) 

./nl 243 ^ gl J gl J v ; 



V. EXAMPLES 



A. The DBI Galileons 



As a first example, we consider the DBI Galileons (2f| described by the action ([5]). Here we extended the original DBI 
Galileon model by introducing V(<p), f(4>) and g{<j>) which depend on </> weakly so that the slow-varying parameters 
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are sufficiently small. 



py, x) = -f(4>r vi - zxJW) + mr 1 - , G ^ x ) = T^kjW) ' (45) 

It is worth mentioning that only particular choices of the functions /(</>) and g(4>) come from genuine higher 
dimensional symmetries, although the equations of motions are still kept to be second order. 

In order to analyse this model, it is convenient to define cd = ^/P,x which corresponds to the sound speed in the 
DBI model in the absence of the Galileon-like term. Since it is known that oc 1 /c 2 D in DBI inflation and we are 

interested in the case where the large non-Gaussianity is generated, we consider only the case with cd *C 1. 

We assume that the inflation is driven by the potential term, that is, V ^> X/cd, gHtfiX in Eq. (fTTj) as in usual 
DBI inflation. During inflation, the field equation (fT2"j) becomes 

±H^ + l^X +v (46) 

CD tfj 

The value of cd is specified from the background equations once V{<j>) and f(cf)) are given. We assume c'd/cdH -C 1 
so that, at the leading order, cd is constant (it is also possible to construct a model where Cd is constant by choosing 
a functional form of f(<j>) for a given V{4>), see Ref. [5l[ for the details). We define a parameter bu = gHcj)/(? D so that 
the first two terms in Eq. (|4"o) are comparable when bjj is of order 1. It is worth noting that bjj can be also expressed 
as bo — v/{e— 2>v). Making use of the concrete forms of P(X, </>) and G(X, <t>) in Eq. (|45l) as well as the definition of 
bo, the coefficients C\, C2, C3 and C4 in the third order action p0|) are given by 



Ci = -V(l + 24M, C 2 = --V(1 + 126 2J ) I C 3 = ~^-b D , C 4 = ^b D , (47) 

where we have used XP XX ~ l/(2c 3 D ), X 2 P XXX ~ 3/(4^), G, x ~ g/c%, XG, XX ~ 2g/c% X 2 G, XXX ~ 6g/c s D 
and cd <C 1. 

Then, the parameters fV> in Eq. (f39|) become 

/( i )= 3 (l + 20b g ) /(2)= (l + 126gj /(3) _ (1 + l2b D )b D (48) 



2cl (1 + 46d)(1 + 126 d ) ' J 4(1 + 46 D )<5, 2(l + 46 D ) 2 c 



D 



From Eq. (|44[) we obtain 

/NL — nntJl ft , /i L \fi , 101. \2 ' V^ y J 



f e<?„« = 5 (21 + 546bu + 3776^ + 60486%,) 
,XL 324 C 2 (1 + 46 D )(1 + 126 D ) 2 



which becomes — 0.16/c 2 for &d ^> 1 and — 0.32/c 2 for bu — > 0. To obtain Eq. (|4T))) . we have used that cd and c s are 
related as 



(1 + 12&B, 2 



(1 + 46 D ) 



<5 • (50) 



The nonlinear parameter f^™ 1 scales as oc 1/c 2 and can be detectable by future experiments such as PLANCK for 
sufficiently small c s . 

Especially, in the case of fj = s = 0, we find that the following relation holds: 

^irl^ 1 -^- (51) 

Combined this with Eqs. (|5"Tj) and (HHJ), we can express f^" a in terms of n s and r as 

fZ U * -20 (1 ~ 2 ns)2 , (52) 

where the coefficient is almost independent of 6d for this set up. This relation suggests that these Galileon-like terms 
do not help embed the DBI inflation into string theory (FH, [53[ . 
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B. G-inflation 

The second example is recently proposed G-inflation where an exact de Sitter solution is realised without introducing 
a potential (42J. In this model the functions P(X) and G(X) are chosen as 

P iX) = -X + ^, G{X) = ±X. (53) 

The de Sitter solution is obtained when Px + 3Hcb/M 3 = is satisfied in Eq. (|TT|) . The solutions are obtained as 

X = l iM 3 x 1 H 2 = — (1 ~ X)2 , where 1 = = = y/E-^- (54) 
^ 18|i i xy/l-x/2 M pl 1 1 

Here as in Ref. [Hj], we only consider simple cases where /j, <C M p i. In this situation, (1 — x) ~ y/S/i/Mpi <C 1 and /it 
is related with M pi , M, i? as fi = 6M^H 2 /M 3 . 

Then, making use of the concrete forms of P(X,<p) and G{X,<f>) in Eq. (j5"3")l . the fine-tuning condition P,x + 
3H(j>/M 3 = and Eq. the primordial power spectrum and tensor to scalar ratio are given by 

where we neglected the correction to the tensor to scalar ratio, r = 16c s e, that arises if we allow a small deviation 
from the pure de Sitter inflation. As was pointed by Ref. [42[ , r becomes nonzero even in the pure de Sitter solution, 
i.e. for e = 0. On the other hand, n s — 1 becomes for the pure de Sitter inflation, but again if we allow a small 
deviation for it parameterised by the slow- varying parameters, n s is given by Eq. (|26[) . 
For the primordial bispectrum, the parameters in Eq. (|39[) become 

2 2(1 — x) (1 — x) 

where we have set x — 1 unless it appears in the form of 1 — x. Then from Eq. (1441) we obtain 

'sr'-ejrbr (57) 

Especially, from Eq. (|55[) and assuming the pure de Sitter inflation, e = 0, f^ q ™ 1 is related with r as 

4 - 62 i» (58) 



r 



which gives a strong constraint on this model. For example if r = 0.17, which can be detected by the PLANCK satellite 
01 j /nl"*' = 15.1. Therefore, a detection or non-detection of the tensor mode and equilateral type non-Gaussianity 
by PLANCK will tightly constrain the model. 



VI. CONCLUSION 

The DBI inflation model [l| has been extensively studied recently for both theoretical and phenomenological reasons. 
Especially, it is well known that the DBI inflation model can give large non-Gaussianity of the Cosmic Microwave 
Background (CMB) temperature fluctuations. Recently, de Rham and Tolley [25[ proposed a new model so called the 
DBI Galileons based on a probe brane action in the higher-dimensional space time. Interestingly, this model naturally 
provides a connection between the DBI model and the relativistic generalisation of the Galileon model [231 where the 
equation of motion is at most second order in derivatives due to the Galileon symmetry d^ir — > d^n + Cp. Since the 
DBI inflation is supposed to be driven by the dynamics of the brane in the higher dimensional bulk, it is interesting 
to study the effect of the Galileon-like terms in DBI inflation 

In this paper, motivated by the DBI Galileons, we studied primordial fluctuations generated during inflation de- 
scribed by the action ^ which is obtained by generalising the action of the DBI Galileons. This generalisation is done 
in a similar way to the extension of the DBI inflation to the K-inflation. In order to calculate the statistical quantities 
of £, the curvature perturbation on uniform density hypersurfaces on large scales, at leading order in the slow-varying 
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approximations, we have adopted the simple procedure [44j to first calculate the bispectrum of the fluctuations of 
inflaton <j> in the flat gauge then relate it to that of ( using the delta-N formalism (fT4")l . 

For the linear perturbations, we have confirmed that, owing to the Galileon-like term, the expression of the sound 
speed c s for the scalar perturbations is modified from the usual K-inflation model (Eq. (|f 71) ). We also provided 
general expressions for the power spectrum V^, spectral index n s and tensor to scalar ratio r at leading order in the 
slow-varying approximations. In these expressions, e defined by Eq. (|22p . not e = —H/H 2 , plays an important role 
(Eqs. (l24|) . (|26p and (|28|l). Due to this, the consistency relation between the tensor to scalar ratio and the tensor 
spectrum index is broken if there exists the Galileon-like term. 

We calculated the bispectrum at the leading order in the slow- varying variables. The Galileon-like term gives a new 
shape A3 in addition to the shapes A\ and A2 which arise in the usual K-inflation (Eqs. (|34p . (|55|) and (|36p ). For the 
new shape ^3, we checked the validity of u sing the estimator for the equilateral- type Non-Gaussianity, f^ q ™ 1 , based 
on the shape correlator introduced by Ref. |48[ and showed that the overlap is at about 99.99% level, which justifies 
the use of this estimator to measure the amplitude of the bispectrum even in the presence of the Galileon-like term, 
we obtained the general expression for the amplitude of the bispectrum in Eq. (I44[) . 

For the concrete examples, we have considered two models: one is the DBI Galileons described by the action ([5]). 
The other is G-inflation model proposed by Ref. [42]. For the DBI Galileons, in the small sound speed limit, f^ lL 
is given by Eq. (|49l) and written in terms of the sound speed c s and bu that is related to the amplitude of the 
Galileon-like term. Since it scales as f^ tl cx 1/c 2 , large primordial non-Gaussianities can be obtained when c s is 
much smaller than 1, similar to the usual DBI inflation. It is worth mentioning that for a given c s , the bo- dependence 
of /^j™ i s wea k an d we obtained —0.32/c 2 . < < — 0.16/c^. For G-inflation where an exact de Sitter solution is 

obtained without any potential terms, we found f^" zl and the tensor to scalar ratio was related as f^™ 1 — 4.62r~ 2 / 3 . 
Although a small deviation from the de Sitter solution could give a correction to the tensor to scalar ratio, this relation 
gives a stringent constraint on the model by a detection or non-detection of the equilateral type non-Gaussianity and 
the tensor to scalar ratio. 

In this paper, we considered the cubic order interaction in the Galileon theory and its relativistic generalisation. As 
is shown in |25| , it is possible to add two more higher order interactions which again arise from the probe brane action 
in a five-dimensional spacetime with the Gauss-Bonnet term. It is also possible to generalise these terms in the same 
way as generalising Eq. (O to Eq. It would be interesting to study phenomenology of this generalisation. The 
single field model arises from a probe brane action in the five-dimensional spacetime. If the DBI Galileons have some 
connections to string theory, the DBI Galileons should be naturally a multi-field model as in the DBI inflation where 
the position of the brane in each compact direction is described by a scalar field. The multi-field Galileon model has 
been extensively studied recently [55M60I ] and the relativistic extension of the model has been proposed Q . We leave 
the study of multi-field DBI Galileons, the relativistic generalisation of the multi-field Galileon, for a forthcoming 
paper. 
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